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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. ' 

TECHNICAL NOTE NO. 331. 

THE RESISTANCE TO THE STEADY MOTION OF 
SMALL SPHERES IN FLUIDS . * 
By R. A. Castleman- 

Introduction. 

There seems to be little reliable information conveniently 
available as to the resistance encountered by small spheres 
moving steadily at moderate speeds in fluids. The present pa- 
per, while presenting nothing new in the way of either theory 

or data, has three objects: first, to show that published data 

i 

are sufficient to furnish approximate information; second, to 
present this information in form convenient for computation; 
and, third, to indicate where further research is needed. 



* This paper was prepared in. the spring of 1924, but the pres- 
sure of other work prevented its being put into shape for publi- 
cation at that time. Meanwhile the work of Liebster and Schil- 
ler, which covers the range of values of ^EL from 0.12 to 

2000, has appeared. Their results, agreeing well with those 
quoted- here, lend confidence to the conclusions drawn in this 

paper, and also close up the gap 200 < ^B. < 3000. Therefore, 

it seems worth while to present this paper, especially since it 
covers a greater range than that covered by Liebeter and Schil- 
ler, thereby allowing conclusions to be drawn that were not ev- 
ident from their work; and it shows in a rather striking way 
how data, apparently unrelated either to each other or to the 
problem under consideration, may sometimes be utilized with prof- 
it by applying to them the principle of dynamical similarity. 
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Theory 

In this paper, the following symbols will be used: 

Table I. 



Quantity 




Symbol 


Dimensions 


Resistance 




X 

*■ 


KLT~ S 


Diameter of sphere 




d 


L 


Head pressure 

Relative speed - 
Sphere vs. L'edium 


D-h 


- P . r 
1 ,S 

±1T d 


ML" 1 T~* 
LT _1 


Density of medium 




P 


ML" 3 


Density of sphere 




a 


ML -3 " 


Viscosity of medium 




P- 


ML" 1 T- 1 


Mass of , sphere 


* 


2.1 


M 


Length of molecular 
free path 


mean 


\ 


L 


Acceleration due to 


gravity 


g 


LT~ 3 



Hydrodynamical theory indicates that, if d and V are 
sufficiently small, and if the fluid may be regarded as (a) 
homogeneous and. (b) infinite in extent, then the flow of fluid 
round the sphere should be laminar, so that the resistance 
should be mostly due to the shearing of the medium r In such a 
case the viscosity |J will be the dominating factor. In fact, 
for these conditions, the law deduced theoretically by Stokes 
(Reference l) seems to have been amply verified by experiment. 
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This lav; is 



3 TT [JL d 7 



(1) 



If, on the other hand, d and V are sufficiently large, 
the resistance should he mostly due to the energy dissipated in 
the formation of eddies "behind the sphere, and hence depend on ■ 
the kinetic energy of the motion. The law of resistance should 
then become: 



where k is a constant. 

The limits of validity of these two simple laws are not 
well known, and there is an intermediate range for which no 
simple mathematical statement of the law of resistance is pos- 
sible. 

lie will see if we can not infer some of this information 
from published data. 

The easiest and most reliable way to compare results ob- 
tained by different investigators, expressed in different units, 
is by the dimensional method. The treatment of this method is 
without the scope of this paper (Reference 2), but our present 
requirements are quite simple- 

If we assume that, with the limitations already mentioned, 
3>k depends only on d, Yy p, and p., then the law of resistance 
must be of the form 



D 



k p d' 



2 Y 2 



(a) 



f d, 7,P,n) = 0 



(3) 
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Dimensional reasoning shows that it mast be possible to put (3) 
in the form . 



_2h__ F (SS£S . JV.-"^(^? V i 3a > 
P v s V ii J 



where the form of the function F must be determined by other 
considerations. For example, in the region where Stokes' law 
(Equation l) is valid, the second member of (3a) becomes . 
In the regime of the square law (Equation 2) it becomes 



Dimensional reasoning, however, tells us more- Its main 
result is the so-called "principle of dynamical similarity." 
In the present case this principle merely says that the value 
Of *H- should be uniquely dependent on that of ^jt&, - that, 
if we form the product in any way whatever, we should get 



the same value of 



Ob 



To investigate the law of resistance, then, we should de- 

termino the form of the function F- To this end, some of the 

published experimental data have been examined and are shown in 

Figs. 1-4, where values of are plotted on logarithmic 

Pt^ 2 

bases against the corresponding values of ME • 

Experimental Data 

I. Silvey (Reference 3) observed the rate of fall of crops 
of mercury, of diameters from 0-012 cm to 0.07 cm (0.005 in. to 



d.YP 
V- 

0.000034 to 0.0066. His results agree with Stokes' law. 



0.03 in.) in castor oil, over a range of values of Q^f- from 
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II. Arnold (Reference 4) observed the rate of fall of Rose 
metal spheres, of diameters fToir. 0-015 cm to 0,14 cm (0.005 in.' 
to 0.055 in.), in colza oil in a tube of internal . diameter 1.09 
cm (0.43 in.). The values of ranged from 0.002 to 2.4. 
Since the proximity of the walls of the fall tube caused a 
violation of condition (b), Page 2 above, the observed speed of 
fall must be corrected to that which would prevail iri an infi- 
nite fluid. Such a correction has been proposed by Landenburg 
(Reference 5) . It is 

*inf. = "obs. (l + 3.4 4) (4) 

where d 1 is the diameter of the fall tube. 

Arnold's observations, corrected by Equation (4), agree 
with Stokes 1 law for the smaller values of Not much 

weight is given to his observations at larger values of 422. 
since in this range his spheres were so large compared to the 
size of the fall tube that eddies could not form freely behind 
the sphere. Hence Stokes 1 law would tend to hold longer than 
it would in an infinite fluid. 

The results are shown in Figs. 2 and 3, where each plotted 

» 

point represents the mean of several observations, the observed 
data being too numerous for convenient individual representation. 

III. Allen (Reference 6) observed the motion of air bub- 
bles and of solid spheres, in water and in aniline^, over the 
ranges indicated in Table II. 
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Table II. 
Observations of H. S. Allen 



Sphere 


Medium 


D iamet er 


dj£ 

M> 


cm 


in. 


Air bubble 
Air bubble 
Paraffin 
Amber 
Steel 


Aniline 

Water 

Aniline 

Water 

Water 


0.00? - 0.11 
0.01 - 0.06 
0.07 - 0.30 
0.12 - 0.30 
0.30 - 0.80 


0.0028 - 0.043 
0.004 - 0.024 
0.028 - 0.118 
0.047 - 0.118 
0.118 - 0.315 


0.009 - 11.0 
0.270 - 25.0 
; 0.500 - 20.0 
21.000 -204.0 
2300 - 8200 



For small values' of for which Ke- used air bubbles 

in aniline. Allen 1 s results are seen to be quite scattered and 
to fallj in" general, below the theoretical curve. The scatter- 
ing is to be expected from the difficulty in measuring the di- 
ameter of the small bubbles (which were caught under a glass . 
and measured microscopically); while the low values obtained 



3 in this range can be accounted for, in part at least, 



for -^4 

by the fact that the air in the bubble was gradually frittered 
away by the medium, - an effect observed by Arnold (Reference 
4). Allen* s measurements of the diameters are thus too small, 
and if corrected for this effect, the points would fall nearer 
the theoretical curve. A.llen f s observations at small values 
of C an therefore not be regarded as contradictory to 

Stokes' law. 



COPY 

IV. Liebster andi, Scn^iiLe'^'-C Reference 7). observed the 

Me to,* oftim u ngle 

Memorial -Aeronautical 
laboratory 
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rate of fall of steel spheres, of diameters 0.1 cm to 0.7 cm 
(0.04 in. to 0.28 in.), in glycerin, sugar solutions, and water 
Their observations covered a range of values of Q$p- from 0.12 
to 2000. Their results are shown in Figs. 3 and 4. 

V» Millikan and his collaborators have investigated the 
motion of very small spheres of various materials, - oil, mer- 
cury, shellac, etc., - in a number of different gases, at val- 
ues of less than about 0.0005- It was found that, as Ion 
as the ratio of the mean free path of the molecules of the gas 
to the diameter of the sphere is very small, Stokes' law holds 
quite well. When, however, this ratio becomes considerable the 
sphere falls faster than is indicated by Stokes 1 law. This oc- 
curs when A. is large (at reduced gas pressure, for instance), 
or d is very small, or both, so that condition (a), Pasje 2 
above, is violated. Equation (l) then becomes 

D = 5 tt n d T (5) 
_ , „ A. 

1+.9 Z 

The results are summarized and the complete form of the 
function q> ^ is given in a paper by Millikan (Reference 8) . 
Under most of the conditions with which the present paper is 
concerned, Millikan' s correction is quite negligible. Its mag- 
nitude is roughly indicated in Table III. 
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Table III. 

Correction to Stokes* Law - Spheres in Air 



Pressure 


Temp. 


Diam. of sphere 


cp . 4 


cm Hg 


in. Hg 


C C 


°F 


cm 


■ in. 


76 


29.9 


20 


68 


0.001 


0.0004 


1.7 


76 


29.9 


20 


68 


0.0001 


0.00004 


16.0 


38 


14.95 


20 


68 


0.001 


0.0004 


3.3 


76 


29.9 


0 


32 


0.001 


0.0004 


1.4 


38 


14.95 


0 


32 


0.0001 


0.00004 


30.0 



Remarks 



The data shown in Figs. 1-4 were obtained by four different 
investigators using three different kinds of material for the 
spheres, - solid, liquid, and gaseous. While only liquids were 
used as resisting media in the experiments there represented, 
Millikan' s experiments verify Stokes' law for motion in gases 
at small values of except under the conditions noted 

above- The author has not found data which seem reliable for 
resistance in gases in the range 0.0005 < < 6000, but 
there seems no reason why this should differ from the curve in 
Figs. 1-4. Hence, while future research should aim to obtain 
information as to the resistance in gases in the range indicated 
above, the following conclusions seem justified! 

1. The quantities involved in Equation (3), with the limi- 
tations indicated on Page 3 above, are the only physical quanti- 
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ties seriously affecting the resistance to the motion of spheres 
in fluids in the range 0.00001 < < 6000. 

2. Stokes 1 law holds fairly accurately for £21 < 0.5, 
and approximately (within about 7$) to ^j^- = 1*0. This means, 
for the latter case, a drop of water of diameter about 0.008 cm 
(0.003 in.) falling in air at 20°0 (68°F) and 76 cm (29.9 in.) 
Hg- 

3. For the range 0.5 < ^jjjS- < 6000, the law of resistance 
is given by the curve of Figs. 3 and 4, with a possible error 

of about 7$. 

There are several facts to which it seems hardly necessary 
to call attention: first, Figs. 1-4 really represent one con- 
tinuous curve, it being broken up into seven parts in the man- 
ner shown for convenience in plotting and computing; second, 
since and ^B- are dimensionless, any self- consistent set 

of units can be at once applied to the relation there shown; 
and, third, the value of the resistance for any particular value 
of i s PV 8 is finally determined, not by the value of d, V, P> 
or u- separately, nor by that of any incomplete combination of 
these, - such as dV, or ^f, - but only by that of the com- 
plete product. 

Applicat ion 
We may illustrate the application of the information con- 



.N.A.C.A. Technical Note No. 231 



10 



taincd in Figs. 1-4, by computing the terminal speeds of spheres, 
of various sizes, falling under the influence of gravity. These 
are obtained by equating the resistance D, to the effective 
weight of the sphere, tt d 3 (o* - P) g. 

Case I (3gB < 0.5} 

i7Td 3 (<7-P)g=3TTU.dV 

Case II (0.5< 3££ < 6000) 

In this region we must use the parameters obtained from Figs. 
3 and 4. To do this, we proceed as follows J 

i.e. __s - gg . i n 43 8 ( ° - p) 



ilT d 2 P V S \ TT d S p V 2 -| TT d S P V S 



_ 2 dg (q - P) = (? v 

- 3 p V s 7 ^ 

Solving Equations' (?) and' (7a), for V and d, we get 



d = SJL (9) 
VP ^ ; 
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Taking simultaneous values of x and y from Figs. 3 and 
4, and substituting proper values of M-,P,o", and g, we my 
compute corresponding values of d and V. 

From Equations (6), (8) and (9), have "been computed the 
terminal .speeds of spheres of unit specific gravity, falling in 
air at 20°C (68°F) and 76 cm (29.9 in.) Hg. The results, ex- 
pressed in c.g.s. units, are shown in Fig. 5. 

Abstract 

Data on the resistance to the steady notion of small 
spheres in fluids, obtained by various investigators, are col- 
lected and presented in logarithmic graphs, the diraensionless 

variables 2h_ and *33E2. being used as coordinates. 
PV 2 V- 

The entire range of values of ^P. from 0.00001 to 6000, 
is found to be satisfactorily represented by a single continuous 
curve. 

Stokes' law holds accurately to =0.5 and approxi- 

mately (within 7$) to = i.o. 

The data are applied to the computation of the terminal 
speed of spheres of unit specific gravity, falling in air. 
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